(S' 

NASA 

Technical 

Paper 

2473 

July 1 985 

Hertzian Contact in Two 
and Three Dimensions 

John H. Tripp 






NASA 

Technical 

Paper 

2473 

1985 

Hertzian Contact in Two 
and Three Dimensions 


John H. Tripp 

Lewis Research Center 
Cleveland, Ohio 


NASA 

National Aeronautics 
and Space Administration 


Scientific and Technical 
Information Branch 



Summary 

The basic solution to the problem of mechanical 
contact between elastically deforming solids was 
proposed by Hertz over a century ago and has been used 
by tribologists and others ever since in a steadily 
increasing number of applications. While the theoretical 
development is not conceptually difficult and treatments 
exist to suit ail tastes, it is nonetheless interesting to trace 
the relationships among the solutions in different 
dimensions. Such an approach is used herein to shed light 
on the curious and sometimes perplexing behavior of line 
contacts. A final object of this brief review is simply to 
gather a number of the more frequently used contact 
expressions together as a convenient reference and for 
comparative purposes. 

Introduction 

In a recent review (ref. 1) Johnson first summarizes the 
fundamentals of contact mechanics as originally worked 
out by Hertz and then considers many of the subsequent 
developments and applications of the Hertz solution. For 
the great majority of cases, an application of Hertz 
theory requires little more than selecting the correct 
formula, and this end could be served by a simple 
handbook compilation. Subtleties do sometimes arise 
when using such expressions, however, particularly for 
line contacts where stresses and strains due to loading of 
a boundary element diminish at large distances only as 
fast as the inverse of the distance from that element. This 
implies that the material displacements themselves 
actually diverge logarithmically with distance, and if no 
boundaries intervene, the result may be hard to interpret. 
For such situations, an appreciation of the derivation of 
the expressions and the underlying assumptions becomes 
valuable. 

Beginning with a fundamental solution to the problem 
of a point force on an elastic half-space, we first show 
how it is used to yield the solution for the problem of 
frictionless contact between the half-space and a rigid 
ellipsoidal punch. This fundamental point force solution 
is then transformed to the solution for a line force acting 
on the half-space, the equivalent of the point force 
problem in two dimensions. The problem of contact by a 
rigid circular punch, a cylinder in three dimensions, is 


then solved in an analogous manner. From this, an 
approximate solution for a roller squeezed between 
parallel flats is derived and compared to the exact result. 
While the conventional Hertz analysis correctly yields the 
contact width for such configurations, its prediction for 
the compliance of nonplanar bodies is poor. This may be 
improved substantially by using elastic solutions 
specifically accounting for the geometry of the deforming 
body. 

Finally, a comparative tabulation of the more useful 
formula in three and two dimensions is given. Since it is 
not one of our present aims, no exhaustive list of 
applications of Hertz theory is provided herein. 
Reference, however, may usefully be made to the 
bibliography compiled in reference 1. 

Symbols 

A 0 semiaxis of ellipsoid in ^-direction, m 

a contact semiminor axis in three dimensions, m; 

contact semiwidth in two dimensions, m 

#0 semiminor axis of ellipsoid cross section in 
contact plane, m 

B complete elliptic integral 

Bq semiaxis of ellipsoid in >’-direction, m 
b contact semimajor axis in three dimensions, m 

b 0 semimajor axis of ellipsoid cross section in 
contact plane, m 

C bounding curve of contact region S 

Cq semiaxis of ellipsoid in ^-direction, m 

D complete elliptic integral; diameter of body in 
^-direction, m 

d compliance, m 

d 0 compliance datum, m 

E Young's modulus, Nm -2 

8 complete elliptic integral of second kind 

F total load in three dimensions, N 

J general complete elliptic integral 

K complete elliptic integral of first kind 

k argument of complete elliptic integral 

p contact pressure in z-direction, Nm~2 
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p mean contact pressure, Nm~2 

Pq maximum (central) contact pressure, Nm-2 

q total load per unit length in two dimensions, 

Nm- 1 

R radial polar coordinate, m; reciprocal of sum of 

principal curvatures of ellipsoid, 
(R ~ 1 + R ~ ] ) “ l , m; radius of circular 
cylinder, m; reciprocal of sum of curvatures 
for two cylinders, (/?f 1 + R 2 ') Hi 
Rq radius of sphere, m 

r radial spherical coordinate, m; radial polar 

coordinate in 5, m 

5 surface of contact region 

u,v,w elastic displacement vector components, m 
x,y,z Cartesian coordinates, m 
z axial cylindrical coordinate, m 

a ratio of principal curvatures of ellipsoid, R y /R x 

7 Euler’s constant 

T ratio of difference to sum of principal 

curvatures of ellipsoid 

e eccentricity of contact ellipse 

f nondimensional z-coordinate, z/a 

6 azimuthal cylindrical coordinate 

k ratio of major to minor axis of contact ellipse, 

b/a (ellipticity) 

v Poisson’s ratio 

£ nondimensional x-coordinate, x/a 

p radial cylindrical coordinate, m 

a normal stress, Nm-2 

r shear stress, Nm-2 

< p azimuthal polar coordinate in 5 

\p azimuthal polar coordinate; digamma function 

Fundamental Point Force Solution 

Any convenient starting point may be adopted in 
deriving the stress tensor to solve a given elasticity 
problem. If the final result satisfies all the general and 
specific conditions of the problem, then the solution is 
acceptable. Thus, the stress tensor for the problem of a 
concentrated point force acting along the normal to the 
plane undeformed surface of an elastic half-space is 
chosen here as starting point, even though it too may be 
constructed from still more elementary results. 

Since the force provides a symmetry axis for the basic 
problem, the coordinates most suited to its description 



(a) 



(a) Cylindrical polar coordinates used for a point force F acting in 
z-direction at the origin on the bounding surface of an elastic half- 
space. 

(b) Plane polar coordinates used for a line force q acting in z -direction 
in x = 0 plane on bounding surface of elastic half-space. 

Figure 1 .—Coordinate systems used in elasticity problems. 

are the cylindrical system (p,0,z ) shown in figure 1(a). In 
reducing the dimensionality of the problem to discuss line 
contact parallel to the y-axis, a system of plane polar 
coordinates (R>\p) in the x,z- plane is convenient, and 
these appear in figure 1(b). Once symmetry is violated, as 
in the ellipsoidal contact problem, the Cartesian system 
(x,y,z) provides a set of fixed directions for 
superimposing stresses arising from parallel but displaced 
point forces. It is straightforward to transform between 
the two polar systems by using the Cartesian system as 
intermediary. 

Consider a point force F acting along the positive 
z-axis on the boundary surface z = 0 of an elastic half- 
space defined by z>0 and having elastic constants E 
(Young’s modulus) and v (Poisson’s ratio). By symmetry, 
polar angle 6 does not appear in the stress tensor and 
components r pl? and t$ z vanish. The four remaining 
components are given by reference 2: 


3 F z3 
° z ~ 2^7$ 

F „ , J 1 z z\ 

-’-T^'-^KrV^ p2r* rl) 

3 F pz 2 
7 ^“ 2tt r5 



where 


r = (x 2 +y 2 + z 2 )i/2 = (p 2 + z 2 ) 1/2 


is the spherical radial coordinate giving the distance from 
the point of application of the force. These stress 
components satisfy the general requirements of 
mechanical equilibrium and compatibility. They diminish 
to zero as the inverse square of r at large distance. On the 
plane boundary, components a z and t pz vanish except at 
the origin, and the net force acting over any 
hemispherical surface centered on the origin is readily 
shown to be Fin the positive z-direction, independent of 
its radius. These are the special requirements of the 
singular solution (eq. (1)) which ensures that it describes 
the correct physical problem. The displacements in the 
directions of increasing (p,0,z) compatible with equations 
(1) are found to be 


Contact In Three Dimensions 

Rigid Ellipsoid and Elastic Half-Space 

Clearly, the point force in equation (1) may be replaced 
by the force p bS acting over a sufficiently small surface 
element bS due to any given distribution of normal 
pressure p on the boundary, and by linear superposition 
the net stresses and displacements will be obtained by 
suitably integrating the resulting expressions. In the 
contact of a rigid punch of any given shape with the half- 
space the problem is to find what pressure distribution 
acting within the actual contact area S will deform the 
plane such that its shape exactly matches the punch 
throughout 5. This matching condition, it turns out, is 
sufficient to determine uniquely the bounding curve C of 
S as well as the pressure distribution within it; hence, the 
stress and strain may be found for any point in the half- 
space. In the absence of friction, the shear stress at the 
surface is everywhere zero. 

To simplify the problem for illustrative purposes, we 
consider matching only the normal displacement 
component w, assuming that frictionless slip 
accommodates the radial movement u within the contact 
region required by the solution. The problem of matching 
radial displacements has been considered by Goodman 
(ref. 3). In the contact plane z = 0, equation (2c) becomes 



pbS 

r 


(3) 


(\-2v)(\ + v)F\z 1 P 2 Z~ 

2t rFp L r (1 -2*) r3 ^ 


(2a) 


v = 0 


(2b) 


f ro +^ 2 

w ~ 2 ^ e [ r3 


+ 


2(1 -y2) 
r 


(2c) 


Clearly, u and w are singular at the origin while the 
arbitrary reference value zero at infinity is approached as 
r~ K To be physically acceptable, the singularities may be 
imagined to be excluded by replacing the point force with 
an appropriate distribution of stress over a finite surface 
enclosing a small volume around the origin. Strictly this 
distribution should be that given by equations (1), 
although, in practice, St. Venant’s principle allows us to 
relax this requirement and to use any distribution which 
yields the same resultant force from the small excluded 
volume. 


by replacing Fwith p bS and introducing the plane strain 
modulus E f =E/( l-^ 2 ). The net displacement at any 
point Min the contact plane can be found most easily by 
setting up the local polar coordinates (r,0) shown in 
figure 2, since in this system bS/r = b<j> br is independent 
of the coordinates. Hence, we have 


Wm== v F i 1 p(r> ^ dr d<t> (4) 

s 

In the neighborhood of the contact, occurring at an 
extremum of the punch with respect to the contact plane, 
the shape may be written by suitably orienting axes in the 
general principal axis form 


2 R x 2 R y 


(5) 
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and the matching condition is simply 



Figure 2. — Local polar coordinates in contact region of loaded elastic 
half-space showing contact bounding curve C. 


where R x and R y are the two principal radii of curvature 
of the punch at the extremum and d is the penetration 
depth of this tip, lying on the z-axis. Later we will further 
assume R y >R x so that the elliptical cross section has its 
major axis in the ^-direction. 

The projection of C on the contact plane is now 
assumed to be elliptical (C in general will not be a plane 
curve) with semiaxes (a,b) along (x,y) while the pressure 
distribution is semiellipsoidal: 

/ x 2 y 2 \ 1 /2 

p<*.y)=Po{ < 6 > 

Integrating p within C yields the total load F : 

F= (7) 


WAf(x>y)=Z p (x,y) (9) 

Performing the r integration first, the integrand may be 
cast in the standard form [a - 20r- yr 2 ] w2 to yield 


l3 ~^[a-2l3r-yr2] l/1 + 
2y 


ay + 02 
2 7 3/2 


arc sin 


(3 + yr 

(ay + 02)1/2 


The 0 dependence of the coefficients (a,0, 7 ) is such that 
incrementing 0 by tt leaves a and y invariant but changes 
the sign of 0. Thus, at the lower limit r- 0, the coefficient 
ensures that the subsequent 0 integration of both terms in 
the interval 0 to will yield zero. At the upper limit of r, 
where the bracket vanishes, the arc sin gives rr/2 and we 
are left with 


WM(x,y) 


Po [ 
4£" J 


2-k 


0 


ay + 02 
7 3/2 


c/0 


(10) 


Examining next the (x,y) dependence of this remaining 
0 integrand, we find a term independent of ( x,y ) and 
terms quadratic in x and y, all of which are invariant 
when adding 7r/2 to 0. A fourth term, bilinear in x and y t 
changes sign. Taking advantage of this property and 
substituting for a and 0 give the result 


» M (x.y)=$ j (^7: 


x 2 sin2 0 
a 2 !? 2 7 3/2 


y 2 

a 2 b 2 


cos 2 0 ^ 

~yFn ) 


d<t> 


( 11 ) 


so that the average pressure F/irab is just two-thirds of its 
central maximum value p$. The consistent determination 
of (a,b) in terms of (R Xf R v ) constitutes a complete 
solution of this problem, since the stress tensor is then 
known in terms of various integrals within C. The 
solution also yields the penetration depth d which in the 
three-dimensional problem may be physically visualized 
as the depth of the tip of the indenter relative to the 
undeformed boundary of the half-space at large 
distances. 

If M in equation (4) is the point (x,y), 


where 7 = (cos 2 0)/a2 + (sin 2 0)/&2. This has precisely the 
functional form required by the matching condition (eq. 
(9)). Thus, the choices for Cand p(x t y) in equation (6) are 
self-consistent, and it at once becomes possible to 
determine the axes of the contact ellipse in terms of the 
principal radii of the indenter. To obtain these in 
standard form, we first introduce two notations: 

K = b/a (12) 

^ \j 0 J2) sin2 0J 3 “7 1/2 03) 


Wa f(x,y) = 


Po_ 
7 r£' 


fi[-( 


x + r cos 0 


( y + r sin 0 V 
\ b ) 


y 

1/2 

dr d(j> 


It is always possible to satisfy k> 1 by appropriately 
naming the (x,y) axes relative to the principal directions 
of the punch. 

Comparing the terms in equations (5) and (11) shows 
the following relationships must hold: 


4 



_J_ = Po 1 ( 

2R X E' Kb\ 
Po 1 


tt/2 


sin 2 (f> 
^3“ 


d<t> 


= Po K [ 

£' 6 J 


ir/2 


COS^(j) 


d<t> ( 14 ) 


» tt/ 2 f X// 2 ( 

cos 2 0 , Po j_ sin 2 j> 

2R y E'nb 3 Q A3 0 E’Kb\ A 


i/<p(15) 


where the second form has been obtained by simple 
partial integration to yield the standard complete elliptic 
integrals B and D , as displayed in appendix A. The 
argument (1 — *-2)1/2 of these functions is just the 
eccentricity e of the contact ellipse. Correspondingly, an 
expression for the penetration depth d is found from the 
term independent of ( x,y ) in equation (1 1): 


d = 


Po & 

E‘ K 




( 16 ) 


In this form d is clearly proportional to the elliptic 
integral K(t). 

To complete the solution it is necessary to determine 
the eccentricity e, or equivalently k, from the known 
principal curvatures of the indenter. In terms of the ratio 
a = R y /R x , equations (14) and (15) may at once be solved 
for k to give 


K = 



(17) 


Substituting for the elliptic integrals B and D in favor of 
K and S A gives 


K = 


[<-■>! 



( 18 ) 


Alternative forms of equations (17) and (18) often 
appear in which the dimensionless ratio a. is replaced by T 
defined as the ratio of the difference to the sum of the 
principal curvatures. Thus, if R ~ 1 = Rx 1 +Ry 1 IS 
sum of the curvatures, then T = (R^ 1 - R ~ [ )R and the 
equations for k become, respectively. 


k — 


( 1+T D\ 1/2 

Vl^r b ) 


and 


k = 


2/w - (1 +T)8] 1/2 

(1-OS . 


(19) 


(20) 


From equation (5), a condition for the contact 
boundary C to be a plane curve is that the radius ratio a 
should be equal to the square of the axial ratio k 2. 
Equation (17) however shows this can hold only if B = D , 
which occurs when e=0, the special case of circular 
contact. In general C is not plane. 

When the punch has penetrated to depth d beneath the 
surface z = 0, its elliptic cross section in this plane has 
semiaxes {a^b 0 ) given according to equation (5) by 


g o _ bp =d 

2 R x 2Ry 


( 21 ) 


From equations (14) to (16) the semiaxes of C satisfy the 
relation 


2R x a+ 2 R y a a a 


( 22 ) 


which in combination with equation (21) leads to the 
interesting and perhaps unanticipated result 


a 2 b2 
«<5 + b 2 0 


( 23 ) 


This simple expression is useful in avoiding the pitfall of 
setting (a,b) = {a 0 ,b 0 ). In the ^ = 0 plane, this distinction 
between a and a 0 is shown graphically in figure 3. 

Adding equations (14) and (15) and substituting the 
total load Ffrom equation (7) yield an explicit expression 
for the semimajor axis of C: 



Figure 3.— ln_y = 0 plane, ellipsoidal punch appears approximately cir- 
cular with radius R x . Its intercept at x = a u with initial contact plane 
Z = 0 is compared with contact semiaxis a. 
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, f 3F«2* !■') /3F,2RS\I'J 

H wr <e+ ° /,,2) J <24 » 

The semiminor axis a = b//c then becomes 



while the penetration depth d is given by 



Since k and R are fixed by the geometry of the rigid 
indenter, independent of load, it follows from equations 
(24) to (26) that the linear dimensions of the contact 
ellipse grow as the one-third power of the load, while the 
penetration depth accordingly increases as the two-thirds 
power. 

For given a, equation (17) or its equivalents (eqs. (18) 
to (20)) are transcendental equations for k which may 
readily be solved, for example, by iteration. When a is 
only slightly greater than unity, e is small and the 
expansions of the elliptic integrals given in appendix A 
may be used to obtain a closed solution. If this is written 
in the form 


K = ct" (27) 


and only the first two terms in the expansions are used, 
the exponent n is easily found to be 2/3. Comparing 
equation (27) with direct numerical evaluation of k shows 
this approximation to be surprisingly accurate, with, for 
example, an error of only 3 percent for a as large as 10. 
For such a radius ratio, e = 0.97 so that the accuracy of 
equation (27) clearly far exceeds that of the 
approximations used for the elliptic integrals. 

In a more exhaustive investigation of the /c(a) 
relationship, Hamrock and Brewe (ref. 4) used linear 
regression to fit equation (27) over the larger range 
1 <a<100. They find the least squares value of n to be 
2/ 7r, only some 5 percent smaller than the limiting value, 
which yields a maximum error in k of less than 4 percent 
for their larger range. 

For the special case of a spherical punch of radius Z? 0 , 
the preceding general analysis is more complicated than 
necessary. The final result, however, is derived readily 
from the general case by setting R=R 0 /2 and T = 0. 
Equation (17) then gives x=l since when argument 
(1 - k-2)1/2 = o integrals Z? and £> both have the value 7r/4 


and K and 8 both equal 7 t/ 2. The radius of the circular 
contact (eq. (24) or (25)) and the penetration depth (eq. 
(26)) become 



These results are the more familiar Hertz expressions, 
written to show that the plane edge of the circular contact 
occurs at a depth of exactly d/2. The actual contact area 
at this level is consequently just half the cross section area 
of the punch at the level of initial contact, z = 0 , which 
also follows directly from equation (23). 

Equations (24) to (26) together with equation (17) 
furnish a complete and rigorous solution of the contact 
problem for a rigid indenter of arbitrary shape on an 
elastic half-space. Substituting the contact pressure (eq. 
(6)) into equations (1) and (2) and integrating over the 
contact region give the stress tensor and corresponding 
elastic displacements at the general internal field point, 
including the plane of initial contact itself. Such 
calculations were first carried out for circular contact by 
Huber (ref. 5) in 1904. The general elliptic contact was 
not analyzed until 1917 by Belajef (ref. 6), almost four 
decades after Hertz’ solution for the surface stresses. 
Somewhat later, Thomas and Hoersch (ref. 7) obtained 
the stress components on the axis of symmetry as a 
function of the eccentricity of the contact ellipse e. Such 
calculations allow the principal and octahedral shear 
stresses to be found; the absolute maxima of these 
stresses determine the initial yield point of the half-space 
according, respectively, to the criteria of Tresca or von 
Mises. The position of the maximum principal shear 
stress, for example, rises from a depth of 0.78a beneath 
the surface when e = 1 (line contact) to 0.57a when e =0 
(circular contact), taking v to be 0.3. The magnitude of 
this stress is, however, insensitive to e, varying only from 
0.30/? 0 in line contact to 0.31/? 0 in circular contact. 

Another feature of interest appears in the radial stress 
component o p which at the surface changes sign from 
compressive at contact center to tensile just inside the 
contact edge. In circular contact this component ranges 
from -0.8/? 0 at the center to 0.133/? 0 at the edge itself. 
For brittle materials, failure is initiated at such points of 
maximum tensile stress. Furthermore, since the 
circumferential stress remains compressive, equal but 
opposite to a p outside the contact circle, the maximum 
shear stress in the surface is quite appreciable. These 
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stresses fall monotonically to zero as the inverse square of 
the distance. 

Extensions of Point Contact 

In most practical applications of this mechanical 
analysis, the axes of the contact ellipse are several orders 
of magnitude smaller than the radii of curvature of the 
indenter; therefore, the solution may be extended to 
cover a wider class of contact problems. Since elastic 
displacements decrease as \/r (see eq. (2)), both bodies 
are then sufficiently flat within the region where 
displacements are significant so that it is immaterial 
whether one or both of the bodies have curvature and 
whether one or both are deformable. The gap between 
the bodies (given by eq. (5)) becomes the sum of the two 
separate deviations from the z = 0 plane, and it is thus 
necessary only to express the gap in principal axis form. 
If the principal axes of the two bodies are aligned, a 
common situation, then, with a suitable convention for 
the signs of radii, this is achieved by writing 
1 /R x = (\/R x ) ] + (\/R x ) 2 and likewise for the 
^-curvatures. Similarly, the elastic displacements of the 
bodies (given by eq. (8)) are summed and \/E' may 
simply be replaced by \/E[ + X/E^ Caution: A 
convention is sometimes used whereby 1 /E' is defined as 
(\/E[ + l/£ > 2)/2, which has certain appeal when both 
bodies are of the same material; however, if this 
convention is adopted, then in the preceding work, E' 
must be replaced by EW 2. Some of the most useful 
results derived in this section for three-dimensional 
contact are displayed in tables 1(a) and (b) at the back of 
the report (pp. 21 and 22). 


Reduction of the Fundamental Solution 

In the previous section the solution for a point force F 
was developed to handle distributed pressure problems by 
replacing Fwith p <5S. A corresponding replacement of F 
by q 6/ allows a treatment of a concentrated line force of 
magnitude q per unit length acting normal to the half- 
space along line element <5/. If q is constant and the entire 
_y-axis acts as the line, an integration of equation (1) 
yields a stress tensor independent of the ^-coordinate. 
Because the p- and ^-directions used in equation (1) are 
referenced to each source element, it is convenient to 
transform the tensor to Cartesian form before carrying 
out the integration. Then, the p- and ^-coordinates of the 
field point M taken in the y-0 plane with respect to 
source F on the y - axis are as shown in figure 4. The 
required tensor transformation, taking account of the 
zero components in the original cylindrical system, 
becomes 



Figure 4.— Coordinates of field point M with respect to point force h 
acting on v-axis. 


o x = a p cos 2 0 + o e sin 2 0 
Gy = o p sin 2 6 + oq cos 2 6 

7 xy = ~ °e) sin 26 

7 XZ ~ 7 pZ cos ^ 

T yz ~ 7 pz ^ 


(30) 


From symmetry it is clear that r xy and r yz vanish upon 
integration, but for the remaining components no way 
has been found to avoid actually carrying out the rather 
lengthy integrations. The results depending only on (x,z) 
may be more clearly expressed in terms of the (R, yp) 
coordinates of figure 1 : 


2q 


o Y = — sin 2 yp cos \p 

7 TR 


2vq 
7 tR 


cos 




7 xy 


= 0 


7 xz = 


2 Q 

7 tR 


sin cos 2 \p 


T .vz = ° 


, z =-^cosH 


(31) 
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It follows that these equations describe a state of plane 
strain since a y = v(o x + <j z ) and r xy = r yz = 0. Note also that 
the only role of the Poisson ratio in this quasi-two- 
dimensional solution is to generate a normal stress in the 
excluded ^-dimension just sufficient to maintain zero y 
strain. The more important involvement of v in the 
original three-dimensional solution of equation (1) is 
required by the compatibility relations. 

Having expressed the Cartesian stress components in 
polar coordinates, it is clear what their cylindrical polar 
counterparts are. By inspection, all components in this 
representation are zero except for a R and a y : 

lq i \ 

o R =-—c°s* 

> (32) 

2v< l ! \ 

Oy = - — cos ^ = va R j 

Thus, in the y-plane, this somewhat lengthy analysis has 
reduced the stress system to a single radial compressive 
stress (with no shear). This system is so simple that it is 
often used as the starting point in elementary treatments. 
The displacements produced by this stress system may be 
found, like the stresses, by integrating the three- 
dimensional result, but it is actually easier to derive them 
directly from the stresses of equation (32). Neither 
approach avoids the problem of the logarithmic 
divergences of the displacements at R~0 and R- oo. The 
former coincides with the singularity in the stress so that 
it makes physical sense to excise a small region around 
the origin in the same manner as for the three- 
dimensional solution. The divergence at large distances in 
regions which are stress-free is somewhat more 
troublesome, since it deprives the solution of an 
undisplaced region of the contact plane to serve as the 
natural reference level for displacements. Instead, a finite 
datum point is forced on the solution and its choice 
affects the form, if not the content, of subsequent 
expressions for various displacements of interest. 

In the cylindrical system of equation (32) one useful 
form of the displacement components ( u,v,w ) in the 
(R,y,\p) directions may be written 

U= ~^E' ( 2cos ^ ln | + l^7^ sin ^) ) 


again showing that this is plane strain. The parameter d 0 
is chosen as the reference depth beneath the surface z = 0 
at which the displacement in the symmetry plane x = 0 
vanishes. The azimuthal component w has been made 
antisymmetric about this plane. 

Contact In Two Dimensions 

Rigid Cylinder and Elastic Half-Space 

For contact problems the important part of the 
solution is again the normal (downward) displacement of 
the original bounding surface z = 0, given in this case by 


w(R, - tt/2) = - w(R, tt/2) = 

2q r, R 1 1 

- .-rL ln 3i + 2ii37, J < 34 > 

This is applied first to the contact of a rigid cylindrical 
punch with an elastic half-space, where initial contact is 
along a generator touching the .y-axis. If only normal 
displacements are matched, the technique follows exactly 
that for frictionless point contact. An ^-dependent 
pressure distribution p(x) dx replaces the line force q and 
the integral of equation (34) is matched to the azimuthal 
cross section of the indenter. Shear stresses in the initial 
contact plane are again zero. The choice for p(x), guided 
by equation (6), is taken to be the semielliptical 
distribution 

p(x)=p (j [\---) (35) 

so that the total load per unit length is just 
Q=\P(A (36) 

The resultant z-component of the displacement in the 
z = 0 plane then becomes 


- w(x) = 


2q_ 

■kE' 




2^0 
7 rE 


>f (•-£) 


1/2 


In I*-*' | dx' 
or, expressing the integral in standard form, 


> (33) 

, . q (- ... R \~2v sin^\ | 

w = -rr, I 2 sin ^ In — ^ cos \p + - — - ) 

7 t £\ d 0 \~v 1 - v / / 
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['" + + 7 )] 


2poa 

ttE' 


d-X’2) 


In x' dx' 

a 


(37) 


When combining equation (34) with the pressure 
distribution of equation (35), r becomes the distance 
between the pressure element and the field point, which 
accounts for the appearance of the modulus of the 
relative coordinate in the integral. Similarly, d is the 
distance between the pressure element and the datum 
point so that equation (34) requires generalization to 
accommodate an arbitrary datum before incorporation 
with the distributed pressure. On the assumption, 
however, that the datum is chosen at a distance large 
compared to the extent of the contact pressure, the 
variation of this distance has been neglected in arriving at 
equation (37). In the three-dimensional case, the 
reference point is taken at infinity and this consideration 
never explicitly arises. 

The remaining integrand in equation (37) is singular 
for points within the contact strip flanking the^-axis out 
to distance a on either side. Taking the Cauchy principal 
value for the integral as indicated by the modulus sign 
leads to the following closed forms: 



As before, the contact dimension must be determined 
by matching the displaced surface to the indenter shape 
within the contact region. Corresponding to equation (5), 
if R x is the radius of curvature of the indenter at the point 
of initial contact, then the local shape in general is given 
by 


z p 



(39) 


For two-dimensional contact, radius R x becomes 
identical with radius R, defined previously as the 
reciprocal of the sum of the x- and v-curvatures. We shall 
therefore suppress the jc subscript in subsequent 
expressions referring to two-dimensional cases. The 
displaced surface according to equations (36) to (38) may 
now be written 


- W(x) = 



In 


26 ^ 

a 


v 

2(T^7) 



(40) 


which matches equation (39) in form. The assumed 
pressure distribution is thus consistent and may be made 
to fit the physical conditions of this problem by choosing 
the contact semiwidth a and penetration depth d to be 



These results correspond to the three-dimensional 
expressions given in equations (24) to (26) and show that 
the linear size of the contact zone increases here as the 
1/2 power of the load per unit length. The penetration 
depth increases essentially linearly with load except, as 
equation (42) shows, there is also a weak logarithmic 
dependence from the appearance of a in the bracket. 
Because of the freedom of choice for the reference depth 
d$, no particular relationship exists between a and d 
which would correspond to equations (26) or (29) for 
point contact. 

Characteristically, do might be assigned a value 
comparable to R. A typical example of the effect of 
choosing different reference positions may be given by 
shifting the datum from (x,z) — (0,tf 0 ), a depth d 0 beneath 
the indenter, to (r/o>0)> a distance to the side in the 
original undisplaced surface. Reexpressing the 
displacements of equation (33) in accordance with this 
new choice leads to +1/2 in place of the term 
- v/2{\ - v) within the bracket of equation (42) for d. 

As the distance \x\ from the contact zone increases, 
equations (37) and (38) show that the normal surface 
displacement increases without limit as In (\x\/d$) in the 
direction opposite to that of the applied force. Such 
behavior is of course expected from the discussion of the 
fundamental displacements given in connection with 
equation (33). 

With the contact width known as a function of load, 
the rigid punch problem can now be solved completely 
for stress and strain tensors by appropriate combination 
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of the pressure distribution of equation (35) with the 
basic line force solution (eqs. (31) and (33)). The first to 
perform such calculations were Huber and Fuchs (ref. 8), 
while more recently, elegant solutions have been provided 
by Poritsky (ref. 9) and Smith and Liu (ref. 10). In 
appendix B, the straightforward integration of equation 
(31) leads directly to the stress tensor given by equation 
(B15). Plots of these stress intensities are presented in 
figures (6) to (8) while some special case:> of interest are 
discussed here in more detail. 

Consider first the surface of the half-space f=0. From 
equation (B12) it follows that a 2 = dt(£ 2 -l) where the 
choice of sign must be determined by the requirement 
that 0 be real in equation (B 1 1 ). This leads to a 2 = 1 - £2 
when £ 2 <1 (inside the contact zone) and a 2 = £ 2 -l 
when £ 2 > 1 (outside the contact zone), so that the actual 
value of a 2 is never negative. Correspondingly, 
0 = (l -£ 2 )i/2 inside and zero outside. When the surface 
is approached by taking lim f— *0, the stress components 
are all found to vanish for points external to the contact 
region. Internal to the contact region, only the shear 
stress t xz is zero, while the normal stresses a x and o z are 
equal to each other and given by 

o x = ° z =-Po( 1-^) 1/2 (43) 

The ^-component, of course, is just the contact 
pressure (eq. (35)), which is one of the boundary 
conditions of the model, but the .v-component is perhaps 
unexpected. It produces a state of two-dimensional 
uniform (hydrostatic) compression in the surface layer, 
as a comparison of figures 6 and 8 also reveals. While the 
vanishing of o z and r xz on the surface outside the contact 
zone are also boundary conditions of the model, the 
vanishing of a x is not a requirement and again contributes 
some new information to the contact problem. In fact, 
this behavior of o x in two dimensions contrasts strongly 
with the behavior of the three-dimensional contact 
discussed earlier, where the corresponding normal 
component of stress in the radial direction actually 
reversed sign and became tensile just inside the contact 
boundary. 

Another special case of interest is found by setting £ = 0 
to obtain the internal stress components on the symmetry 
axis of the indenter. In this case, a = 0 = (1 + f 2 ) 172 . Again 
the shear stress t xz vanishes, as we expect by symmetry. 
The nonzero normal components are given by 

a x = - [2f2+l_2rtl + r 2 ) 1/2 | ] 

( 1 + 72 ( 

> (44) 

a = - — 1 

(l+f2) 1/ 2 ) 


which are plotted in figure 5. The normal stress a v , equal 
to v times the sum of o x and a z> is also shown, again for 
the value of y = 0.3. At large depths these components 
display explicitly the 1/f distance dependence 
characteristic of the plane strain condition. 

Now that the normal stresses are no longer equal, the 
maximum shear stress r m lies in planes bisecting the x-and 
^-directions with a magnitude of \o z -o x /2. This 
component is displayed in figure 5, which shows the 
previously noted maximum value of 0.3 p 0 at a depth 
approximately three-fourths of the Hertz semiwidth a , 
the initial yield point. By contrast, since r w vanishes 
everywhere at the boundary, the material immediately 
adjacent to the indenter is dead and not expected to flow. 

Extensions of Line Contact 

Extension of the analysis to cases where both 
contacting members have curvature and elasticity is less 
trivial in two than in three dimensions, when it sufficed 
simply to introduce a composite curvature and elastic 
modulus for the pair. The validity of this simple 
adaptation rested directly on the effective vanishing of 
the displacement in regions where curvature produced 
significant departure of either undeformed body from the 
ideal plane assumed for the basic solution. In two 
dimensions no part of this surface remains undisplaced so 
that the procedure of transferring curvature from one 
side to the other, known as bending the gap, represents at 
best a crude approximation which might be expected to 
break down in some situations. 

Consider first the contact of an elastic cylinder with a 
rigid plane, the inverse of the situation described by the 
fundamental solution. If E' and v now refer to the 
cylinder, then equation (41) still gives a good 
approximation to the contact width provided as always 
this turns out to be small compared with R. In the 



Figure 5. — Stress components for points along indenter axis internal to 
elastic half-space. Both distance and stress are given in nondimen- 
sional forms defined by equations (Bl) and (B2). 
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expression for d (eq. (42)), however, the depth is now 
measured in the cylindrical body rather than in the half- 
space, and at distances d 0 -R from the contact the shape 
of the cylinder may depart significantly from that of the 
half-space. Moreover, the parabolic form of equation 
(39) is often taken to approximate a circular roller, in 
which case the infinite half-space has been bent into a 
space finite in two dimensions. This emphasizes the 
nature of the approximations involved and the need to 
exercise judgment in selecting a datum. A natural choice 
in such cases is to fix the center of curvature by setting 
d$ = R, whereupon the distance d of equation (42) 
becomes the shortening 8R of the radius in the direction 
of the applied force q\ 




v 

2{T^V)_ 


(45) 


The approximately equals sign in equation (45) is a 
reminder that this estimate of 3R rests on the as yet 
untested bending hypothesis. A more accurate expression 
is presented in equation (47) based on an elastic solution 
for the line force acting on a circular cylindrical surface. 

For parallel contact of two different cylinders, it 
follows from equations (39) and (40) that the contact 
semiwidth a may still be obtained from equation (41), as 
in three dimensions, by interpreting R~ 1 as (tfr'+flf 1 ) 
and E' ~ 1 as (E\ ~ 1 + E^~ '). To find the approach of 
centers, however, these variables alone are insufficient 
because of the appearance of R and v in equation (45) 
referring individually to each cylinder. The compliance is 
most clearly expressed just by adding the individual 3R 
values in which, of course, a and q are common. 
Reduction to a single equivalent roller and flat is, in this 
respect, no longer particularly simple. 

For the diametral reduction of a roller in compression 
between two parallel rigid flats, it is merely necessary to 
double the right side of equation (45). If the flats 
themselves are elastic, however, a little care is required 
first to substitute appropriate composite variables when 
finding a and second to appreciate that if the flats are of 
dissimilar materials, the two contact widths will not be 
equal. Thus, diametrically opposed radii are shortened by 
different amounts whose sum gives the net compliance. 

As a final application, this analysis will be used to 
demonstrate its own limitations in treating the 
deformation of any body other than the half-space. 
Consider again the preceding example of the diametral 
compliance of a cylinder between rigid flats. If the gap 
bending hypothesis is correct, the analysis should yield a 
diametral compliance independent of the arbitrary 
reference point. Let this datum be chosen some distance 
1<R from the midplane along the shortened diameter. 
Then, according to equation (42), the compliance 3D 


obtained by summing the decrements of the two segments 
will be 



(46) 


This expression has explicit, though weak, dependence on 
the datum via /. Yet another result is obtained if d$ = 2R 
is inserted directly in equation (42). Use of the bending 
hypothesis for two-dimensional problems involving 
compliance is thus generally inaccurate or even 
inconsistent (ref. 1 1). 

A resolution of this difficulty is to begin with a 
different basic elastic solution applying directly to curved 
bodies. One such solution, for a circular cylinder loaded 
by two equal and opposite line forces acting at the ends of 
a chord, has been given (ref. 12) and used in a contact 
analysis (ref. 13). This leads again to the pressure 
distribution of equation (35) in which the contact 
semiwidth a of equation (41 ) is, however, reduced slightly 
by the appearance of the factor (1 + Iq/nE'R) in the 
denominator. From the discussion of equation (42) it 
might be expected that the effect of the exact solution on 
d would be seen chiefly in the second term in the bracket, 
already shown to be sensitive to choice of datum. Thus, 
in addition to requiring the new expression for a in the 
first term, the rigorous solution also replaces - v/2{ \ - v) 
by (In 2 - 1/2): 

sr - jT? + ( , n2 -o] ,47> 

Setting a!R- 10 - *, the fractional error in the earlier 
approximation for a amounts to (\0~ n /2) 2 whereas for 
3R given in equation (45) it is far larger, roughly -0.2 /n 
for typical values of v. In practical situations n could be 
as small as 2; thus, in the worst case, application of plane 
Hertz theory to a cylinder overestimates a by a negligible 
0.0025 percent but underestimates 3R by as much as 10 
percent. 

It may be seen then that the gap bending conjecture 
works well in determining contact widths but encounters 
difficulties in predicting compliances where the arbitrary 
datum is involved. From this it is perhaps not surprising 
that the contact semi width of equation (41 ) may be found 
as the appropriate limit of equation (25) for the 
semiminor axis in three-dimensional contact, whereas the 
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corresponding link between compliances in two- and 
three-dimensional situations does not exist. The 
connection between line and point contact loaded to the 
same maximum pressure p 0 may be established by 
eliminating p 0 between equations (7) and (36) to obtain 


Q 


3 F 
2 2b 


(48) 


As R v increases without limit, soT-1, k- oo, 8 — 1, and 
A'-* In 4 k. To maintain a finite load per unit length of the 
elongated contact ellipse, the total load F also diverges 
while satisfying equation (48). Substituting these limiting 
values in equation (25) yields the correct two-dimensional 
result for the contact semiwidth (eq. (41)). However, 
while the corresponding limit of equation (26) correctly 
predicts a logarithmic divergence in the two-dimensional 
compliance, it clearly cannot be used to introduce a 
datum level d 0 . Tables 11(a) and (b) exhibit graphically a 
number of the more important results of this section on 
two-dimensional contact problems (see pp. 23 and 24). 

Line Contact of Real Surfaces 

A further interesting development of two-dimensional 
contact analysis was carried out by McCallion and 
Truong (ref. 14) when they investigated the effect of 
surface roughness on the load-compliance relation for a 
rough cylinder between smooth flats. As a model of the 
pressure produced by the asperity layer they used the 
family of contact pressures given by 

( *2 V 

P=Po (^-^ 2 ) \x\<a ( 49 ) 


acting in opposition over a band at each end of a 
diameter of the circular cylinder. The contact analysis is 
based on the elastic solution for diametrically opposed 
line forces acting on a cylinder (see appendix C) so that it 
is limited to a/R<< 1 from the outset. The standard 
integral corresponding to that of equation (38) now 
appears with the variable exponent 0 in place of 1/2 but, 
again provided a/R<< 1, the diametral compression can 
still be evaluated in closed form. The result equivalent to 
equation (45) is found to be 


bR 


- % I'-- 

wE L Q 


2 * . - 0 
— + In 2 - 


0+1/2 


As with the effect of transferring curvature from the rigid 
to the elastic contact member, the effect of changing the 
distribution of load within the narrow band of width 2a is 
confined to the second term of equation (42). For the 
special case 0 = 1/2 required by the matching condition in 
smooth contact, ^( 1 )= -7 and this second term reduces 
to In 2 - 1/2 in agreement with the result shown earlier 
(eq- (47)). 

At this point a number of results for the compliance of 
a cylinder and flat configuration have been derived, such 
as those found in equations (42), (47), and (50). The 
important differences in these expressions are determined 
by which of the two contacting members possesses 
curvature or elasticity, by choice of the datum, or by 
distribution of contact pressure. Further modification 
may be produced by matching tangential displacements 
as part of the condition of contact. Continuing in this 
vein, a paper by Nikpur and Gohar (ref. 15) exhibits no 
fewer than nine different compliance formulae extracted 
from the literature for comparison with a tenth 
expression based on their own experimental data. Beyond 
a number of purely numerical errors, probably arising 
from notational difficulties, most of the apparent 
conflicts between theories clearly arise from genuine 
differences in the mathematical models. These are 
precisely the type under consideration herein, concerned 
with whether to include compliance of the platens, where 
to fix the datum if this is included, what contact pressure 
to use, and so forth. 

Although all ten expressions predict a compliance 
increasing with applied load, as required for mechanical 
stability, there is no such agreement for the dependence 
on cylinder radius which ranges from monotonic 
increasing to constant and even to decreasing. Despite 
this, the work of reference 14 does illustrate a further 
important point beyond simply the need to recognize 
what assumptions are built into a given model. This 
sometimes overlooked point concerns the physical 
validity of the assumptions themselves under the actual 
conditions of measurement. The behavior of a real 
surface often turns out to be rather different from that of 
its ideal, smooth, isotropic counterpart (ref. 16), 
demanding considerable skill and ingenuity in the 
construction of faithful rather than merely solvable 
mathematical models. It is this fact which continues, 
some 100 years after the Hertz analysis, to justify both 
theoretical and experimental studies of such an 
apparently straightforward mechanical system. 


Conclusions 



\p(j3 + 1/2) + y 


(50) 


One objective of this review has been to examine some 
of the assumptions entering into the analysis of 



mechanical contact between elastic bodies. Particular 
emphasis has been placed on reducing three-dimensional 
solutions to quasi-two-dimensional situations where 
elementary metric properties of space introduce 
difficulties by virtue of the logarithmic divergence of 
elastic displacements with distance from the contact. The 
variety of forms this permits for the compliance of 
contacting bodies in two dimensions underlines the need 
to exhibit unambiguously the assumptions on which each 
rests. 

The application of plane Hertz theory to the 
displacement of cylindrical surfaces is based on the gap 
bending hypothesis which asserts that the effect of 
geometry on mechanical conditions in the contact region 
is determined entirely by the algebraic sum of the 
curvatures of the two contacting bodies. Stated another 
way, contact between two curved surfaces can always be 
reduced to contact between a plane and an equivalent 
curved surface, whose curvature is the sum of the two 
original curvatures. Predictions of plane Hertz theory 
applied to cylinders have been compared to those of an 
analysis based on exact solutions for the displacement of 
cylindrical surfaces by a line force. As with its optical 
counterpart in lens theory, the bending conjecture 
produces results of mixed quality. Thus, in two 
dimensions, the contact width is predicted with high 


accuracy but the hypothesis breaks down for the 
compliance, resulting in errors in practical cases as large 
as 10 percent. If for a particular application only the 
contact width and the semielliptical pressure distribution 
are required, then, as argued by Johnson, the errors 
introduced by the plane Hertz approximation are of 
comparable magnitude with errors incurred by using 
linearized elasticity theory, and the extra effort required 
by the refined solution is hardly justified. In matters 
involving compliance, however, the argument is invalid 
and it becomes important to improve on the plane elastic 
solution, if for no other reason than to provide at least 
some error estimate of the initial approximation. 

While the elastic analysis presented herein should be 
sufficiently general to allow modification appropriate to 
any configuration, it nevertheless seems worthwhile to 
collect some of the more common results together for 
easy reference. Accordingly, tables are included as a 
useful and illustrative summary of the principal results of 
this review. 


Lewis Research Center 

National Aeronautics and Space Administration 
Cleveland, Ohio, August 1984 
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Appendix A 

Complete Elliptic Integrals 


General Definitions 


Series Expansions 


A general complete elliptic integral is defined by the 
expression 


For e very small (nearly circular contact): 


J tt/2 

jd<t> 

o A 

where A = (l —k 2 sin 2 0)1/2. For particular forms of 4>, 
the conventional notation is as follows: 


4 1 3 

-B= 1 + -(2+ — f 4+. 

7T o 64 


4 _ , 3 , 15 . 

-D= 1 + -e 2 + e 4 4- . . . 
7r 8 64 


-8=1 - -e 2 - 


64 


^ 4 — , 


4> 

Notation 

COS 2 0 
sin 2 0 
A 2 
1 

B 

D 

8 (second kind) 
K (first kind) 


Only two of these elliptic integrals are independent since 
it is clear that B + D-K and B + (1 - k 2 )D= 8. The 
argument of the elliptic integrals appearing in equations 
(17) to (20) and (24) to (26) is given by k = e where 
e = (1 - *-2)1/2 is the eccentricity of the contact ellipse. 


2 1 9 

-K=l+ t€ 2 + 7^4+... 
7r 4 64 


For k very large (long, narrow elliptical contact) 

b=i -K a -0"- 2 -- 

£> = (A-1)+ ^ (a-*)k-2 + ... 

8=i+ K a_ 0' c " 2+ - 

K = A + 1 (A — 1)k-2 + ... 


where A = In 4k. 




Appendix B 

Stress Tensor for Indentation of Elastic Half-Space by Rigid Cylindrical Punch 


Equation (31) represents the stress tensor for a 
concentrated normal line load q on the surface of an 
elastic half-space, using the coordinates of figure 1(b). In 
plane strain, two of the shear components vanish while 
the normal stress in the suppressed direction is just v 
times the sum of the other two normal stresses, leaving 
only three components to solve for. Under contact 
conditions, the line force q is replaced by the pressure 
p(x) given by equation (35) acting over the infinitesimal 
width dx y the net result being expressed as the integral 
over the whole contact region where p(x) is nonzero. 

Coordinates (R,\p) are relative to the position of the 
element dx so that the natural system to use is the 
Cartesian system of figure 1(b). Introducing the 
nondimensional coordinates 


— = — 2 


a-r 2 ) 1 / 2 (s-n 




The result may be written 


3/,/a 


t xz \ = *- 


3/,/ar 


and the nondimensional stress tensor 


To evaluate I\ we note it may be cast in the form 


o 

o — — 

Po 




the integrals for the three components of interest become 




(£-n 2 n 


where z is the complex variable £+/f. Written in this 
way, the integral in equation (B8) is now' of standard 
form and may be expressed as 


= - \ J (1 | 4 2)1/2 « -nr 2 | # 1 (B3) 


I o-r 2 ) 1/2 


' = W \ Z ~ 


tf-iy 


\^U)j \ f 3 / 

where now R2 = (£ - £ ' )2 + {-2. 

The integrals involved here may all be evaluated from 
the single integral 


(1 — £ ' 2 ) 1^2 
—w—W 


and its partial derivatives 


Note that / 2 (z) is analytic everywhere except when \z\- 1 
at the edges of the contact region on the bounding 
surface, so that it is actually the z-derivative of 


h(z)=\ (l-{' 2 ) ,/2 In (*-{')<«' 

J - 1 


In turn, / 3 is closely related to the integral in equation 
(38) which essentially is the limit of / 3 as f— 0. Thus it is 
possible to derive all the present two-dimensional results 
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from / 3 (z) as in reference 9, although it proves 
algebraically simpler to begin with l { (z). 

To obtain / j we first extract the imaginary part of 
(z 2 - 1) 1/2 , denoted by 0. The actual expression for 0 may 
be written 


- ^ 



j 3=[ a 2-(^2-^-i)] 1/2 /V2 


(Bll) 



(B14) 


with a given by 


More explicitly in terms of a and 0, equations (B14) may 
be written 


a 4 = (£2 — £2— l)2 + 4£2f2 (B12) 

In terms of the auxiliary variable 0 we have 

/t = *G3/f-l) (B13) 

leading to 


o x = 2f- 0 - ({2 + p + 1 + a 2 )(p/20a 2 ) j 

?xz = (Z 2 + P- 1 - <* 2 )(tt/20a2) > (B 1 5) 

a z = - 0 + (f 2 + f2 + 1 + a 2 )(p/20a2) ] 

In figures 6 to 8 these stress components have been 
plotted on planes at various depths f beneath the contact 
plane as functions of the transverse coordinate £. 



Figure 6.— Curves of stress component o x perpendicular to load q 
applied along normal to surface f = 0. 




Appendix C 

Two Opposed Line Forces Acting on a Circular Cylinder 


The stress tensor for the problem of two diametrically 
opposed concentrated line forces acting on a circular 
cylinder is derived readily from that of equation (32) for a 
single line force on a plane boundary by observing that 
the surfaces of constant radial stress in that case are 
circular cylinders tangent to the half-space at the origin 
with diameter D = R / cos \j/. The corresponding radial 
stress value is - 2q/izD. Superimpose on this the stress 
arising from an opposed line force q } acting at the 
opposite end of some diameter D (see fig. 9). Since 
^ + i/q = 7r/2 for any point M on this cylindrical surface, 
then for such points only the (R,\l/) directions for q 
coincide with the (\p\jR j) directions, respectively, for q\. 
Thus, if q = q j, the state of stress at M in the x,z-plane 
becomes an isotropic pressure (~2q/D). Finally, reduce 
this pressure to zero by superimposing everywhere an 
isotropic tension of the same magnitude. The cylindrical 
surface is now stress free except along the lines of action 
of the two forces themselves. The cylinder may thus be 
mechanically isolated from the elastic half-space, and the 
resulting stress tensor becomes the required solution for a 
cylinder of the particular diameter D. For the general 
interior point with coordinates (R,\l/) or (R\,\fo), the stress 
components referred to fixed (x,y,z) directions have the 
plane strain values 



Figure 9. — In jc,z-plane two concentrated forces act along normal to 
bounding surface of elastic half-space, shown shaded. Force q acts at 
surface along + z-axis while force q\ acts at depth D in opposite direc- 
tion. M is any point on circle constructed with segment D of z- axis as 
diameter, having polar coordinates (R,\p) and relative to 

points of action of two forces. 


2 q( 1 cos j/' sin 2 cos i/^ sin 2 ^ ^ 
° x = V \D R / 

2vq ( 2 cos f cos \ 

a y =v{a * +a ^-V\b--R rTJ 

_ 2q ( 1 cos 2 ^ cos3 \j/ { ^ 

° z ~ vU” r rT ~ / 

2 q ( cos 2 \p sin ^ cos 2 sin ^ ^ 
Txz= ~ V V R i R , / 


>(C1) 


T vz ~ 7 xy ~ ® 


Corresponding Cartesian displacement components are 
given by 


u = - 


JLiz^r, 

7 tE ' 1 - V L 


W + lh)- 


1 

2 ( 1 - 2 .) 


x (sin 2\p 4- sin 2\p^) - ^(R sin \p + R\ sin ^ t ) 


v = 0 


w= - 


Tt 2 In ~ + -- - (sin2 i p - sin2 

7 TE |_ R] \ —V 


1 -2v 1 

p(R cos \p-Ri cos t j) 


i]) 


_Q_ 

ttE 


(^ X „D - dp \-2vD-2d 0 \ 
' \ d 0 1 -v D ) 


(C2) 


In this representation, the datum point still lies on the 
Z-axis at distance d 0 from the application point of the 
first force q. Since the shortening dl(x) of a chord parallel 
to the z-axis at distance x is just the difference between w 
values at its ends, it is clear from the separability of terms 
in do shown in equation (C2) that 6/ is now quite 
independent of reference point. When x<<D , 
approximate values for {R,i>\ R\ t ^\) are (x,7r/2; D,0) at 
the upper end and (Z>,0; x,tt/2) at the lower end of the 
chord. The shortening is found to be 
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(C3) 


by substituting in equation (C2) for w. 

The logarithmic singularity in equation (C3) is easily 
handled in the usual way by distributing the concentrated 


line force across a small but finite width. A distribution 
of the type in equation (49) then leads to the compliance 
result of equation (50). In particular for 0=1/2, the 
contact matching condition is met for compression of the 
cylinder between parallel flats, provided the contact 
semiwidth is obtained from equation (41). The 
compliance is then given by equation (47). 
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TABLE I— THREE-DIMENSIONAL STRAIN 


I Rigid bodies are shown shaded. When one body only is elastic, its elastic modulus is E\ while subscripts are used to distinguish bodies of 
different modulus. A similar convention applies to radii of curvature. Negative values denote the radius of a surface concave as viewed from 
outside the material.] 

(a) One body elastic 


Rigid ellipsoid on elastic half-space a 

General ellipsoid, semiaxes A lh B {h C ( , Sphere, radius R 0 
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Rigid flat on elastic ellipsoid or sphere 


All results for the rigid ellipsoid on elastic half-space apply, d is the approach 
of the center of the elastic bodv toward the flat. 


Corollary: 

The diametral shortening of an ellipsoid or sphere compressed between two 
parallel rigid flats is 2d. 















TABLE L— Concluded 


(b) Two bodies elastic 


Two elastic ellipsoids with aligned principal axes 



R x (« V )| + (/?,); 


1 1 1 

E' ~ e; 

In terms of these composite variables all formulae of table 1(a) apply. 
d = d\ +d 2 is the approach of centers, where d, = [E' /E']d is the compliance of 
body i for i~ 1,2. 


Corollary 1: 

For contact between elastic ellipsoid (1) and flat (2), take the limit as the radii 
of body 2 become infinite, d is the approach of center of ellipsoid 1 toward 
undeformed surface of flat, d\ is the compression of the ellipsoid itself; and di 
is the indentation depth of the flat. 


Corollary 2: 

For compression of ellipsoid between two dissimilar ellipsoids or flats, 
appropriate upper (w) and lower (/) composite variables are required, yielding 
individual upper and lower values of quantities a, b, and d. Net diametral 
compression for the central ellipsoid 1 is [E^d u + E/d/]E{. 


Corollary 3: 

For contact of cylinders crossed at right angles, take the limit as </? v ) 2 and 
(/?„) 1 become infinite. 












TABLE II.— PLANE STRAIN 


(Rigid bodies are shown shaded. When one body only is elastic, its elastic modulus is £' , while subscripts are used to distinguish bodies of 
different modulus. A similar convention applies to radii of curvature. Negative values denote the radius of a surface concave as viewed 
from outside the material.] 


(a) One body elastic 


Rigid cylinder on elastic half-space 
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Rigid flat on elastic circular cylinder 3 



The contact semiwidth is given by equation (A). 

The radial shortening or approach of center of elastic body toward the 
flat becomes 


d = 


2q_ 

x£' 


(-" - 0 
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Corollary: 

The diametral shortening of a cylinder compressed between two 
parallel rigid flats is Id. 


a 1« obiain approximate results tor noncircular cylinders, R must he assigned two distinct meanings: to find contact semiwidth a. R is the radius of curvature at the initial contact 
point; to find compliance d. R is the distance of the geometric center of cylinder from the initial contact point. The accuracy of this approximation for d decreases as these two values for 
R diverge. 
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